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Abstract. We study the constraints aris- 
ing on the expansion parameters c and 
d of the Pion electromagnetic form fac- 
tor from the inclusion of pure space-like 
data and the phase of time- like data along 
with one space-like datum, using as in- 
put the first derivative of the QCD polar- 
ization amplitude 77 '(—Q 2 ). These con- 
straints when combined with other anal- 
yses, provide a valuable check on a de- 
termination of c due to Guo et al. and 
on our previous work where pionic con- 
tribution to the (g — 2) of the muon was 
used as the input. This work further illus- 
trates the power of analyticity techniques 
in form factor analysis. 



1 Introduction 

The Pion form factor continues to be of current interest [H 
ElEHUEUHj . In [2j[6] , we developed a framework for obtain- 
ing constraints on the low-energy expansion coefficients c 
and d of the Pion form factor using data from the space- 
like region (t < 0) [TllMJlllOlIll] , and the simultaneous 
inclusion of the phase of time-like data and one space-like 
datum respectively, where c and d are the Taylor coeffi- 
cients in the low-energy expansion of the Pion EM form 
factor given as, 



F*(t) = 7^(0) + -(rl)t + ct 2 + dt 3 + 
o 



(1) 



In [5] , we used as input the pionic contribution to the (g — 
2) of the muon. This technique can be fruitfully extended 
to other inputs. 

Our work has been greatly inspired by the results of 
Caprini [12j who has studied in great detail the problem 
of obtaining allowed regions in the c-d plane using the 
first derivative II'(q 2 ) of the QCD vacuum polarization 



amplitude 77 (q 2 ), which satisfies the following dispersion 
relation. 



// 2n _ 1 (°° lran(t + ie) 



n'( q 2 ) = - f 

T JO 



(t - q>f 



dt, 



(2) 



first with no constraints from data, followed by the inclu- 
sion of the phase of the form factor up to an energy de- 
noted by ti n . Caprini has shown that the modulus of the 
form factor in the time-like region can further improve the 
bounds, via an extensive construction of a function that 
is analytic in the cut plane, where the cut begins at tj ra . 

In the present work, we wish to apply the methods 
developed in Ref. OH to 77 '(q 2 ). Firstly, it was pointed 
out in Ref. [B], that if c is taken to lie in the range 4.49 ± 
0.28 GeV -4 , a value presented in Ref. [13) based on chi- 
ral perturbation theory (for important early work on the 
Pion electromagnetic form factor in chiral perturbation 
theory see Ref. p3]), a possible conflict arises between the 
regions isolated in our work using the phase of time-like 
data and one space-like datum where (g — 2) is taken as 
the input, and that of Caprini that uses magnitude and 
phase of time- like data and 77 '(q 2 ) as input. Therefore it 
is contingent to check whether such a conflict also arises 
if we were to use II'(q 2 ) itself as the input. One of the 
objectives of this work is to explore this issue. Secondly, 
a very recent determination of c yields the value 4.00 ± 
0.50 GeV -4 [TS]. We must therefore ask what the status 
of the purported discrepancy is for this set of values. It 
turns out that the discrepancy is substantially mitigated 
for values of c in the range (3.5 GeV -4 - 4.0 GeV -4 ). In 
fact with the earlier determination of the same quantity 
3.2 ± 0.5 ± 0.9 GeV -4 [TBJ, the central value is accommo- 
dated in the overlap region of the ellipses we find in the 
present work, and the one isolated by Caprini with mod- 
ulus as well as phase of time-like data. In particular, one 
may conclude from here that the value of d is typically 
about lOGcV -6 for c - 3.2 GeV" 4 and about 20 GeV" 6 
for c ~ 4.0 GeV -4 . This is one of the important conclu- 
sions of this work. 

This paper is organized as follows. In Section[5]we pro- 
vide a general discussion on the origin of bounds on which 
the current analysis is based and present an extensive dis- 
cussion of the formalism. This is followed by a discussion 
on the inclusion of space-like data and our results in Sec- 
tion [3l Section |4] discusses the inclusion of the phase of 
time-like data and one space-like datum and our results. 
We then present our final conclusions and a discussion in 
Section [5l 



2 Formalism 

The Pion electromagnetic form factor (t) enters several 
observable quantities through expressions of the type: 

1 



4mJ 



dtp{t)\Fn{t)Y 



(3) 



where p(t) > in the region of integration. Consider as an 
example the first derivative of the QCD vacuum polariza- 
tion 11(—Q 2 ) evaluated at a space-like value of q 2 = — Q 2 , 
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which satisfies the dispersion relation given in Eq. [2] and 
hence considering only the Pion contribution as in [12] . 
n'(—Q 2 ) satisfies the following inequality: 



I7'(-Q 2 ) > 



1 



48tt 2 



dt 



(t + Q 2 ) 2 



1 - 



t 



3/2 



\Mt)\ 2 , 



(4) 

where t n = 4m \. Since this expression has a positive- 
definite integrand, it would be of interest to find a lower 
bound for it. 

This objective is achieved by first mapping the form 
factor in the t plane to the z plane through the following 
conformal map: 



z + 1 
z-1 



t t-n 



^71 



(5) 



resulting in a unit circle and the branch cut lies at its cir- 
cumference. As a result F^(z) is analytic within the unit 
circle in the z plane. Note that this definition is in accor- 
dance with the convention of Caprini [12] and is related 
to the map defined in J5J[S] through z ^ -z. As a result 
of this transformation, Eq. [?] becomes, 



1 r 2 * 
I > — / dO\h(expit 
2?r J 



(6) 



where 



H z ) = f(z)wA z ) 



(7) 

and z = exp(i0). Here, / is the value of Q 2 ), f(z) is 
the form factor in terms of the conformal variable and 
w 7r (z) is the "outer function" to be discussed further. 
Within the unit circle, f(z) is analytic, therefore, the func- 
tion h(z) admits an expansion given by 



h(z) = a,Q + a\Z + aiz 2 + 



(8) 



where the a n are real and are functions of the expansion 
coefficients c and d of the Pion EM form factor given in 
Eq. [T] From the Parseval theorem of Fourier analysis, the 
integral in Eq. [5] is now simply given by, 



n=0 



(9) 



Due to the fact that this quantity is a sum of squares, 
truncating the series yields a lower bound for I. 

In the past, the observable that was studied in great 
detail was the pionic contribution to the muon anoma- 
lous magnetic moment (g — 2), i.c, I was the bound on 
the (g — 2) of the muon, where the normalization and 
the charge radius of the Pion were supplied from exper- 
iments, and the series was truncated after the first two 
terms [17] . A corresponding outer function for the rele- 
vant dispersion relation gives rise to a set of coefficients 
a n with n = 0, 1, ■ • ■ . One could have improved the bound 
/ by retaining the expression up to the three terms shown 
above, and supplying in addition the second Taylor coef- 
ficient coefficient (c). If we were to retain the series up 



to the first four terms, we would have had to supply the 
third Taylor coefficient (d) as well. Improvements on the 
bounds were made possible by using experimental infor- 
mation on the Pion form factor at space-like [18] and at 
time-like [19] values of the momentum transfer. 

In this work, as in [2][6], we ask the reverse question: 
given the value for /, which in this case is the value of 
II'(—Q 2 ) and the available data in both space-like and 
time-like regions, can one obtain bounds on the Taylor 
coefficients c and d of the Pion form factor? 

In Ref. [5], we considered only the inclusion of pure 
space-like data taking up to three constraints. Our results 
were encouraging and provided constraints that were sig- 
nificantly stronger than those found by Caprini |12j using 
only phase of time-like data and the value II Q 2 ). 

In Ref. [6] , we considered the problem of wiring in the 
constraints from the phase of the form factor in the time- 
like region and the constraint from one space-like datum. 
The corresponding Lagrange multiplier technique was de- 
veloped and applied. One significant finding was that the 
simultaneous inclusion of these constraints isolated the al- 
lowed region as one that was significantly smaller than the 
mere intersection of the constraints taken one at a time. 
In [6], we have also carried out a detailed analysis with 
several alternative models for the phase of the form fac- 
tor, using the simple scattering phase formula, two Roy 
equation fits which, via the Fermi- Watson theorem, get 
related to the phase of the form factor, and also directly 
from the parametrization of the form factors themselves. 

Specializing to II'{— Q 2 ), which obeys the inequality 
in Eq. [4] the map from the i-plane to the z-plane results 
in the following closed expression for the outer function 



w-x{z) = 



16 



1 1 



Qnt^ (1 — zd^ 2 



where, 



cL = 



y/U + Q 2 - Vu 



(10) 



(11) 



Retaining only a finite number of terms N = 3 in Eq.[5] 
the expansion coefficients for the function h(z) in terms of 
f{z) and Wtt(z) may be readily expressed as follows, using 
the expansion for the Pion form factor (Eq. [T]) , 



Ol 



ao=h(0)=w*(0), 



h'(0) = <(0)- -rlt„w„(0), 



(12) 



(13) 



a 2 = 



h"{0) 



2! 



«>*(o) ( 



32 cti 



2<(0) ( --rlU 



<(0) 



(14) 
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and 



tion, 



h"'(0) 1 

= = g MO) (-12r^ ff + 384c4-384d4)] 



3! 



+ 



3<(0) 



-r 2 t w + 32ct 2 



+ -[-2<(0)rft r + <'(0)]. 



(15) 



Scaling Eq. ^ by / and keeping only N terms, we get 



A' 



Mo 



E( c «) 2 < i. 



n=0 



where, 



(16) 



(17) 



Any c„ which satisfies /^q < 1 is allowed and the equality 
gives the bound. We have already seen that this yields an 
ellipse in the c — d plane. 

The bound I is evaluated at Q 2 = 2 GeV 2 so that 
I = 0.009546 GeV" 2 and r 2 = 0.42 fm 2 , as in Ref. [T2], 
in order to have a meaningful comparison. Including up 
to the second (third) derivative for F^{t) results in con- 
straints for c (c and d). 



3 Space-like constraints 



Additional constraints in the c — d plane are now obtained 
by wiring in information coming from the space-like re- 
gion, which can be expressed as linear constraints: 



h(x) 



0, 



(18) 



where h(x) includes the value of the form factor F„(t) at 
space-like points defined by Eq.[7]and x is real, lying in the 
range (0 < x < 1) as a result of the conformal map. The 
constraints are included through the method of Lagrange 
multipliers. We set up the following Lagrangian: 



- Vc 2 
2 ^ 

n=0 



^a™^-^^), (19) 



m— 1 



n=0 



where J m = h(x m )/Vl, h(x m ) are the space-like con- 



straints defined at 



1,2,3- ■-. In Eq. {19 



we wish to consider only finite number of expansion coeffi- 
cients c„ and therefore we set N = 3. With the space-like 
constraints included, the equation satisfied by the coeffi- 
cients c„ (analogous to Eq. I16[) is obtained by eliminating 
the Lagrange multipliers a m and can be written conve- 
niently in the form of the following dctcrminantal cqua- 
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1 X<x x\ x\ (1 — X2X1)' 



(1-X1X2) 



= 0. 



(20) 

Note that we use the scaled coefficients c n as opposed to 
a n in Ref. [2] for the sake of notational convenience and as 
a result we need to scale h(x m ). Solving for c and d gives 
an ellipse in the c — d plane, when / is used as an input. 

The data we use are given in Tables [IJ [2] and [3] that 
list the values of x(t), the space- like points for a given t 
and the corresponding value of F n (t). The value of h(x) 
can be evaluated using Eqs. ITlandfTOl To make our nota- 
tions clear, we refer to the data points corresponding to a 
particular |f| as x\, X2 and so on, in the ascending order 
of magnitude of \t\. The tables also show the experimental 
errors in the data. 

The data from various experiments are chosen so as 
to give reliable numerical results. As in Ref. [2] we have 
carried out a detailed study of the bounds obtained by 
considering up to three space-like constraints, both one 
at a time and simultaneously. As expected, the larger the 
number of constraints, the smaller is the allowed region in 
the c — d plane. For one space-like constraint, data com- 
ing from the lower \t\ region give more stringent bounds 
compared to data from the higher \t\ region. On the other 
hand, data from higher \t\ allow us to include a larger 
number of constraints. While this sensitivity is true for 
raw data from experiments, if we were to obtain theo- 
retical fits for the data as an intermediate step it would 
be possible to use low \t\ values as well as include more 
constraints. 

We now present the results of our analysis that are 
captured in several figures. Fig.[T]shows the allowed region 
in the c — d plane obtained by implementing one space- 
like constraint from the Amcndolia data, corresponding to 
the point x\. It may be seen that the range for c is now 
considerably reduced to (—55 GeV -4 , 40 GeV -4 ), whereas 
having no additional space-like constraints yields a range 
of (-210 GeV -4 , 200 GeV -4 ) for c [12]. The corresponding 
ranges with (g — 2) with and without additional space-like 
constraints arc (-20 GeV -4 , 15 GeV -4 ) and (-45 GeV -4 , 
85 GeV -4 ) respectively, where we see similar narrowing 
down of the range for c. One immediate conclusion is that 
the space-like data alone are not expected to constrain c 
significantly for either of the inputs. Therefore, the simul- 
taneous inclusion of space- like constraints and the phase of 
time-like data, which is the subject of the next section, be- 
comes interesting. In Fig. [21 the result from the constraint 
with X2 is presented. As this corresponds to a larger value 
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Table 1. Spacelike data from Bebek et.al [9] 



t(-Q 2 ) [GeV 2 ] F*(t) x(t) 

1 -0.620 0.453 ± 0.014 0.499 

2 -1.216 0.292 ± 0.026 0.606 

3 -1.712 0.246 ± 0.017 0.655 



Table 2. Spacelike data from Tadevosyan et al. [TT| 



t(-Q 2 ) [GeV 2 ] F„(t) x(t) 

1 -0.600 0.433 ± 0.017 0.494 

2 -1.000 0.312 ± 0.016 0.576 

3 -1.600 0.233 ± 0.014 0.645 



Table 3. Spacelike data from Amendolia et.al [10] 



t(-Q 2 ) [GeV 2 ] 




x(t) 


1 -0.131 


0.807 ± 0.015 


0.242 


2 -0.163 


0.750 ± 0.016 


0.275 
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Fig. 1. The allowed ellipse with Amendolia data corresponding 
to Xl. 
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Fig. 2. The allowed ellipse with Amendolia data corresponding 
to X2- 
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Fig. 3. The allowed ellipse with Amendolia data correspond- 
ing to xi, when the data is varied over the experimental un- 
certainty. 



of \t\, the allowed region is enlarged, as expected. In Fig. [3] 
we present the results obtained with one space-like point, 
when the datum corresponding to x\ is varied over its ex- 
perimental errors. Some sensitivity is seen, as the value 
of \t\ is relatively small. In contrast, it may be seen from 
Fig. H] that if a larger value of \t\ is chosen, the sensitivity 
to the experimental error is lower. 

In Ref. [2] we had presented a very detailed discussion 
on the bounds arising from space-like constraints alone, as 
these were found to be stringent when (g — 2) of the muon 
was used as the observable. Since the corresponding con- 
straints with II ' (q 2 ) as the input are not as stringent, it is 
not very fruitful to investigate them in such great detail. 
Nevertheless, we show in Fig. [S] the results arising from 
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Fig. 4. The allowed ellipse with Bebek data corresponding to 
xi, when the data is varied over the experimental uncertainty. 
Note that the dependence on \t\ is rather weak as Bebek data 
lies in the high \t\ region. 



implementing three spacelike constraints simultaneously 
using the data set of Bebek et al. It isolates a region that 
is at variance with other determinations. In particular the 
value of c favoured by chiral perturbation theory appears 
to be excluded. However, since we are not carrying out a 
detailed error analysis, this result need not be considered 
very significant. The use of three simultaneous constraints 
from the Tadevosyan data leads to a result where there is 
no such discrepancy (see Fig. [6]). As a next step, we explore 
the simultaneous inclusion of the phase of time-like data 
with one space-like datum in the following section. The 
results can be easily extended to include more space-like 
points, which is beyond the scope of this work. 



4 Phase of time-like data and one space-like 
datum 



> 

u 
C 



-3000 




Fig. 5. The allowed ellipse with Bebek data corresponding to 
1 (xi), 2 (xi,X2) and 3 (xi, X2, X3) constraints. As the num- 
ber of space-like constraints are increased, the ellipse shrinks 
indicating tighter bounds for c and d. 



> 




Fig. 6. The allowed ellipse with Tadevosyan data correspond- 
ing to 1 (xi), 2 (x\,X2) and 3 (xi, X2, X3) constraints. The 
largest one corresponds to the one space-like constraint and in- 
termediate ellipse to two space-like constraints and the smallest 
to three space-like constraints. 



Time-like phase is the argument of the form factor along 
a part of the cut. Consider the phase of the form factor 
in a region t w < t < t- m in the complex t plane. Assuming 
that this coincides with the phase of the elastic two-pion 
scattering phase shift jT2] we have, in accordance with the 
Watson final-state theorem, 



Avg[F n (t + ie)]=Sl{t), t n <t<t b 



(21) 



The conformal map, as defined in Eq. takes t„ — > z = 
e l9 = -1 9 = 7T and t in -> z = e lB = z in => 6 = 6- m . 
The upper and the lower edges of the cut in the t-plane is 
mapped onto the lower and upper unit disc in the z-planc 
respectively As a result, 



Arg [F„(re l 
Arg [F^re 1 



Sl(9), 9 
51(9), n< 



1 < 9 < 7T 

< (2tt - L 



0, 



(22) 
(23) 



where 5\(9) = 5\(t(9)) and t(9) = t n + Ucotg 2 (9) , which 
is the conformal map with z = cxp(i(9). 

The phase is introduced through the Omncs function, 
denoted as O^(z) in the z-plane, by, 



On (z) = exp 



2jt 



d6 



sl(e) 



1 



ze l 



(24) 



Since the phase of the form factor along the cut (i.e. < 
t < ^in) is compensated by the phase of the two pion 
scattering phase shifts, the following condition holds, 



Im lim 



1 



F^re 1 



= 0, 



(25) 



which is the relevant time- like phase constraint. Substi- 
tuting the expansion for F w (z) in terms of h(z) and the 
outer function Wn(z) the constraint equation becomes: 



]T c„Im lim [[W(9)Y 1 r n e m9 \ = 0, 



(26) 



where, W(9 = C = e lB ) = w n (9)O n (9) 
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Consider now the simultaneous inclusion of the above 
constraint (Eq. |26|) and the constraint from one space-like 
datum: 



J{z) - C nZ n = 0. 



(27) 



where, J(z) = h(z)/\/l is the space-like datum mapped 
on to the z plane, as defined in Eq. I is the bound 
from 77 '(-Q 2 ). 

The space-like and time-like constraints are included 
through the method of Lagrange multipliers where we set 
up the following Lagrangian 

n=0 

oo 

+ -Vc„lim / X(9)\W{9)\Im\\W(9)}- 1 r n e me }d9 

7T '-^ r-»l r 
oo 

+ a(J(z)-Y / c n z n ), (28) 

71 = 

and eliminate the unknown multipliers a and X(9) from 
the Lagrange's equations and obtain an expression for the 
bound. 

The equation for X(9) was derived in Ref. [6] and is 
presented below for completeness, where we assume that 
the first N coefficients are already constrained through 
the normalization of the form factor F^(0) and the Pion 
charge radius r„. We therefore have: 



N 



= -X(9) + J2- 



n=0 



sm(n0-$(0)) 



1 



+ 



d9'X(9 



1 sin [(TV + l/2){6 - 9') - <P{9) + $(9')] 



sin 



2 J 

1-z 2 



+ - d9'X(9 / )(l-z 2 )(3(9)(3(0') + J^ Tr P(9), (29) 



where 



(3(6) 



sin [(TV + 1)6> - <P(9)} - z sin [N6 - ${6)] 



1 + z 2 - 2z cos(6>) 
The expression for a is 



(30) 



2\N+1 



N 



J-J2 C nZ n 



Z 



N+l 



d9'X{9')f3(9') 



We get for the following expression for the bound: 

N , JV 



(31) 



2 
ft) 



Y,(c n ) 2 + ^Y,c n f d9X{e)sm[n6-${9)} 

N \ 

J-J2 c nZ n <1- (32) 



n=0 



In our case, we will set N = 3 following [T2], where r = 
9 : 6 in < 9 < 2ir - 9 m . We solve for A(0) using Eq. ([29]), 




/ 



time-like phase 



time-like phase 
+ 1 space-like 



-100 -50 50 
c [GeV 4 ] 



-50 
c [GeV 



50 



-100 -50 50 
c [GeV 4 ] 



Fig. 7. Bounds on the expansion coefficients c and d. Left 
panel shows the bounds obtained from one pure space-like da- 
tum, the middle panel from the pure phase of time-like data 
and the right from the phase of time-like and one space-like 
datum taken together. Note that the simultaneous inclusion 
substantially shrinks the allowed ellipse. 



obtain the corresponding value of a from Eq. (|3 1 1) . Using 
Cyj's, already constrained by normalization of F n and pion 
charge radius, ^ is evaluated (Eq. (J32)). Only those coef- 
ficients, c„, which satisfy /1q < 1 are retained. Our general 
formalism has been cross-checked against related results 
presented in [20] for the pi K system with N = 2. 

The space-like datum we use for our results are the 
data points corresponding to t(—Q 2 ) = —0.6 GeV 2 from 
the Tadevosyan data set (Table [2]) and t(-Q 2 ) = -0.131 
GeV 2 from the Amendolia data set (Table [3]). They are 
referred to as "Tadevosyan" and "Amendolia" respectively 
in the text and figures. 

The time-like phase is defined as [T2"] . 



8\ (t) = arc tan 



and 



where m„ 



r P (t) 



9Wf 



m 2 — t 



4m; 



(33) 



t 



3/2 



(34) 



770 MeV is the mass of the p meson, T a 



150 MeV is the width of the p resonance and = 139 MeV 
is the mass of the Pion. At low energies Eq. ([3"3"| agrees well 
with the one-loop chiral perturbation theory expression 
for the two-pion elastic scattering phase shifts and also 
with experiments for t > 0.5 GeV 2 , as noted in [TJ] (see 
also Ref. [21]). Therefore, we assume that the phase of the 
pion form factor coincides with Eq. ([3"3"]1 for t n < t < t ln , 
where t- ln = 0.8 GeV 2 . We also use the results for the phase 
shift coming from the Roy equation analysis of Ref. |22] 
given in Eqn. (D.l) therein, for two choices of input pa- 
~) = (0.225,-0.0371) (ACGL) and 



[ao,a 



rameters that are I 

(0.220, -0.0444) [23] (CGL). For most of the work, we have 
chosen U n to be 0.8 GeV 2 in order to have a meaningful 
comparison with the work of Caprini. 

In Fig. [7] we present the bounds when the phase of the 
time-like data with the analytical model for the phase shift 



Gauhar Abbas, B. Ananthanarayan, S.Ramanan: Pion EM form factor and i7'(-Q 2 ) 



7 




" 25 
> 



1 1 1 1 1 1 1 1 1 1 1 

- Min 


' 1 1 1 1 1 1 1 ' 
Central f 


1 

Max •• 


-1,1,1,1,1,1- 


-1,1,1,1,1,1- 


/ Tadevosyan 

-1,1,1,1,1,1- 



- + + 

-2 02468 -2 02468 -2 02468 

c [GeV 4 1 c rGeV 4 1 c [GeV 4 1 



Fig. 9. The allowed ellipse with Tadevosyan datum and phase 
of time-like data, when the datum is varied over its experimen- 
tal uncertainty 
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Fig. 10. The allowed ellipse with Tadevosyan datum and phase 
of time-like data from ACGL [22] 



(Eq. and the space-like datum from Tadevosyan. The 
weak bounds obtained from the space-like data (left panel) 
is already improved by the inclusion of time-like phase 
alone (central panel). When the data from the space-like 
and the time-like region are used together, the allowed 
ellipse shrinks significantly as seen in the right panel in 
Fig. [71 Note that the pure time- like phase constraint ellipse 
shown in this figure was already presented in Ref. [12] . 
while the bounds from pure space-like datum is shown in 
Fig. [5] In Fig. [5] we present a close up view of the allowed 
region for t- m = 0.8 GeV 2 as well as for t- m = 0.6 GeV 2 . 
We observe that the bounds do exhibit some sensitivity 
to the choice of <; n . In order to test the sensitivity of 
the results to other inputs, we now vary the space-like 
datum over its allowed experimental range. The result is 
presented in Fig. [5] In order to judge the dependence on 
the parametrization of the phase shift, we use the two Roy 
equation fits ACGL and CGL [22[[23] . The results for these 
are shown in Figs. [TOl and [TT1 From Figs . [9l [TOl and [TTl we 
can conclude that the bounds are sensitive to the errors 
in the space-like data, while they are not significantly af- 
fected by the parametrization for the phase shifts. 

In order to study the dependence on \t\ we consider the 
datum from the Amendolia data. The bounds (not shown 
here) are very well constrained, however, such small \t\ 
data lead to results that may be highly sensitive to the 
errors in the experiment and consequently are unreliable 
unless we use a theoretical fit for the data, which is beyond 
the scope of our current work. Hence we do not display 
these results. 

It may be concluded from this analysis, that the si- 
multaneous inclusion of the phase of time-like data and 
one space-like datum leads to values of c in the range (—2 
GeV -4 , 8 GeV -4 ) while the range for d is in the range 
(-75 GeV" 6 , 80 GcV~ 6 ). In particular from Fig. [HI it is 
observed that the latter is significantly more constrained 
when II'{q 2 ) is used as the input instead of (g — 2) of 
the muon. The ranges of c in both the cases have a con- 
siderable overlap, which safely accommodates the various 
prior determinations of these quantities in chiral pertur- 
bation theory It is worth emphasizing that the interplay 
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Fig. 12. Comparing the results for the bounds on c and d 
using a M (left panel) [6] and 77 '(g 2 ) (right panel) as inputs. The 
space- like datum comes from the set of Tadevosyan in Table 
while the time-like phase uses the model given in Eq. [33] 



of the form factor. In [SJ, we pointed out that the phe- 
nomenological fits to the form factors from ALEPH also 
yield a number for d ~ 10 GeV" 6 that agrees with our 
determination at lower c values. Another discussion may 
be found in [H], where (c,d) = (3.9 ± 0.1 GeV" 4 , 9.7 ± 
0.4 GeV" 6 ) are reported. 

Our results for the bounds on the Taylor coefficients 
using time-like phase as well as one space-like datum can 
be extended to include the modulus of the form factor 
along a part of the cut. Combined error analysis of space- 
like and time-like data will yield predictions for c and d 
with reliable error estimates, which could be considered in 
the future. 

BA thanks Department of Science and Technology, Govern- 
ment of India for support. SR thanks the Abdus Salam Inter- 
national Centre for Theoretical Physics, Trieste, Italy, for its 
hospitality when part of this work was done. 



of space-like and phase of time-like constraints together 
have narrowed down c and d to an extent that one would 
not have anticipated from the results of either taken sep- 
arately for n ' . 



5 Discussion and Summary 

In this paper, we study the improvements on the bounds 
of the low-energy Taylor expansion coefficients of the Pion 
EM form factor, when both phase of time-like data and 
one space-like datum are used. We use the method of La- 
grange multipliers to include the constraints. In contrast 
to Rcf. [2][6] we have used 77'(— Q 2 ) as an input. This was 
used by Caprini [12) . but without the use of space-like 
constraints. 

We have found that many of the earlier determina- 
tions [2lf6|ll2] are reinforced by the present work, in par- 
ticular the determination of the allowed region for c. How- 
ever, each of these studies provides a different allowed re- 
gion in the c — d plane, although there are regions of in- 
tersection. As in [BJ, we can try and inspect the region of 
intersection of the ellipse in e.g. Fig. [8] with that of Fig. 3 
of [12] that uses the phase as well as the modulus of the 
time-like data along a part of the cut. It appears that for 
c - 4.5 GeV -4 , there is no longer an overlap. However, 
for smaller values such as c~ 4.0 GeV -4 and ever more so 
for c = 3.5 GeV" there is substantial overlap. However, 
a complete picture can be arrived at only when the lat- 
ter analysis is subjected to tests on the sensitivity of the 
results to variations in 7;„ as well as model dependence. 

We find it heartening that so many different approaches 
have led to a coherent picture for these important expan- 
sion coefficients of the pion electromagnetic form factor. In 
particular, the value of d is typically about lOGcV" 6 for 
c - 3.2 GeV" 4 and about 20 GeV" 6 for c - 4.0 GeV" 4 . 
The value of d for c ~ 3.2 GeV" 4 agrees with the that 
from [12] using constraints from phase as well as modulus 
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